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PREFACE 
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1. INTRODUCTION 

The 1978 Agricultural Soil Moisture Experiment (ASME) conducted at a site near 
Colby, Kansas, was designed to acquire field data for use in testing and 
evaluation of soil moisture profile models. Weekly measurements of soil water 
content as a function of depth were made by a combination of .gravimetric 
sampling and neutron-meter techniques. Four locations in each of fourteen 
16.2-hectare (40-acre) fields were monitored from May 20 to the end of 
August 1978. The soil type for all fields is reported to be Keith silt loam ’ - 
(refs. 1, 2). 

Most of the soil moisture models selected for testing and evaluation require 
data on soil hydrologic properties. 1 A detailed knowledge of these properties 
is essential for the analysis of soil water movement. Although field-measured 
soil hydrologic properties are emphasized in soil moisture studies, laboratory 
data on undisturbed core samples are relatively easier to obtain. In some 
instances (e.g., Rogowski, ref. 3; Nielsen et al., ref. 4), predictive models 
utilizing limited measurements and soil structural properties have been found 
to give satisfactory results. 

Colby soils~data that were collected on the physical properties include: bulk 

density, saturated hydraulic conductivity, and a few data points for the rela- 
tionship between soil water pressure and soil water content measured on 
disturbed (crushed and sieved) soil samples. Since hydrologic properties 
determined on disturbed soil samples are not considered suitable and saturated 
hydraulic conductivity represents| only one data point for the relationship 
between hydraulic conductivity and soil water content, information available 
on Colby soils was considered inadequate. Consequently, predictive procedures 
were adopted to develop the needed information. In predicting the hydrologic 


^he relationships between (a) soil water pressure and soil water content and 
(b) hydraulic conductivity and soil water pressure or soil water content. 
These relationships are typically described by curves which span over the 
range of wetness of interest. 
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properties of Keith silt loam, use was made of the available data (refs. 1,2) 
on some of the basic characteristics of the soil. 

This report presents the hydrologic properties of Keith silt loam soil 
obtained from predictive models. Characteristics of the models which were 
used are discussed, and input requirements and computation procedures are 
outlined. The relative merits of the regression and analytical approaches are 
commented upon. Additional efforts for developing and testing predictive 
procedures are recommended. 
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2. AVAILABLE DATA 
2.1 AGRICULTURAL TECHNOLOGY COMPANY 

Bulk density, soil water pressure versus soil water content relationship, and 
saturated hydraulic conductivity data for Keith silt loam collected by the 
Agricultural Technology Company, McCook, Nebraska, are shown respectively in 
tables 2-1, 2-2, and 2-3. 

Bulk density data were single values for the various depths in each of the 
14 fields and showed wide variations between fields and depths (see fig. 2-1). 
Since bulk densities are known to vary over short spaces, more representative 
values for each field were obtained by pooling the data on the basis of 
similarity. Assuming all fields were in Keith silt loam, it was argued that 
differences in surface bulk densities would be more a result of historical 
cultural treatments (e.g., tillage and cropping) and that differences in the 
lower layers would reflect the effects of spatial variabilities in the factors 
of soil formation. On this basis, bulk densities from the surface down to a 
depth of 38 centimeters (15 inches) for all fields of one type — corn, 
pasture, wheat, or fallow — were pooled together, forming a separate group for 
each type of field. For each group, a single smooth curve was hand drawn to 
represent the bulk density profile for the 38 centimeters (15 inches) of 
surface soil. For the lower layers, fields were grouped together on the basis 
of similarity in the pattern of bulk density changes with depth, and smooth 
curves were drawn through the pooled data. 

Field bulk density data were limited to a depth of 137 centimeters (54 inches), 
whereas the water content measurements extended down to a depth of 183 centi- 
meters (72 inches). This necessitated extrapolation of bulk density profiles 
below the 137-centimeter (54-inch) depth. Bulk density data interpolated from 
smooth curves are shown in table 2-4. 

Data for soil water pressure versus soil water content relationships 
(table 2-2) were obtained on disturbed (crushed and sieved) soil samples; 
hence, they were not considered entirely satisfactory. Numerous studies 
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(e.g., Croney and Coleman, ref.. 5; Sharma and Uehara, ref. 6) have shown that 
water retention in the wet range is strongly influenced by soil structure. In 
the dry range, the influence of soil texture is more dominant. Since crushing 
destroys the soil structure, data obtained on such samples cannot approximate 
field conditions. They are useful, however, in the range of soil dryness 
where soil structure has least influence. The range of wetness over which the 
influence of soil structure is significant is somewhat difficult to define, 
for it varies from soil to soil and is modified by the sample history. 
Typically, however, soil water pressures from 0 down to -3000 centimeters can 
be considered to represent the wet range for most soils. Thus, only Colby 

data for lower pressures could be considered useful. 

■ . . » 

Saturated hydraulic conductivities (table 2-3) were obtained at various times 
during a 48-hour flow. The data show fluctuating conductivity values, but 
such fluctuations are not uncommon. Incomplete saturation and structural 
changes during the prolonged wetting and flow process can have both negative 
and positive influences on the flow rate. It is doubtful, therefore, whether 
a mean of all observations for the 48-hour period would be a more reasonable 
value than that of the first few observations. 

Since the objective of the experiment is to closely approximate the field 
structural condition, it is generally recommended (e.g., Klute, ref. 7) that, 
after complete saturation, conductivity should be obtained when approximately 
100 milliliters of water has passed. In the case of Colby data, the length of 
saturation time and the amount of water collected for each observation are not 
known. However, in view of the probable structural disturbance that could 
have taken place, a mean of the first three observations (a 4-hour flow) was 
considered appropriate. 
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2.2 SOIL CONSERVATION SERVICE SOIL SURVEY LABORATORY 

In addition to data supplied by the Agricultural Technology Company, descrip- 
tions of two Keith silt loam profiles (ref. 8) were available from the Soil 
Conservation Service (SCS) Soil Survey Laboratory, Lincoln, Nebraska. * These 
profile surveys were made in 1958. From these surveys, soil texture, organic 
carbon, bulk density, and 15-atmosphere water content data were extracted and 
are reported in table 2-5. For the two profiles, variations in soil texture 
with depth appeared very similar; therefore, particle-size data were combined 
to reflect the variations with respect to soil horizons only. Particle-size 
distributions for the A, B, and C horizons are shown in figure 2-2. The data 
show that the C-horizon [i.e., soil below the 102-centimeter (40-inch) depth] 
is distinctly coarser than the A and B horizons. Although the textural 
compositions of A and B horizons appear somewhat similar, owing to variations 
in structure and organic matter content, water retention characteristics of 
the two horizons can be significantly different. Therefore, treating the two 
horizons separately would be appropriate. 

As expected, organic carbon decreased with depth. These data are plotted in 
figure 2-3, and a smooth curve is hand drawn. 


*An agency of the U.S. Department of Agriculture (USDA). 
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TABLE 2-1.- BULK DENSITY DATA FOR KEITH SILT LOAM 
[As reported by the Agricultural Technology Company] 


Field 

Depth, 
. in. 

Bulk 

density, 

gem" 

Field 

Depth, 

in. 

Bulk 

density, 

.'3 

gem J 

Field 

Depth, 

in. 

Bulk 

density, 
gem" 3 

1 

3 

1 

6 

3 


11 

3 

1 


3* 

»i!« 


3* 



3* 

mmmm 


5* . 



5* 



5* 

MEM 


8 

1.17 


8 



8 

1.13 


9* 

1.34 


9* 

1.40 


9* 

1.43 


15* 

1.45 


.15* 

1.43 


15* 

1.46 


28 

1.22 


24 

1.10 


28 

1.31 


54 

1.09 


52 

1.51 


52 

1.11 

2 

3 

1 .26 

7 

i 

3 

1.39 

12 

3 

1.12 


3* 

1.15 


3* 

1.25 


3* 

1.27 


5* 

1.17 


5* 

1.29 


5* 

1.42 


8 

1.22 


8 

1.25 


8 

1.03 


9* 

1.23 


9* 

1.36 


9* 

1.40 


15* 

1.43 


15* 

1.35 


15* 

1.37 ‘ 


28 

1.40 


28 

1.27 


28 

1.47 


54 

1.27 


52 

1.29 


52 

1.39 

3 

3 

1.34 

8 

' ! 3 

0.94 

13 

3 

1.29 


3* 

1.20 


3* 

1.17 


3* 

1.12 


5* 

1.29 


5* 

1.15 


5* 

1.21 


8 

1.22 


8 

1.14 


8 

1.10 


9* 

1.36 


9* 

1.27 


9* 

1 .34 


15* 

1.38 


15* 

1.34 


15* 

1.36 


28 

1.59 


26 

1.51 


28 

1.34 


54 

1.25 


54 

1 .47 


52 

1.23 

4 

3 

1.09 

9 

3 

1.39 

14 

3 

1.06 


3* 

1.18 


3* 

1.25 


3* 

1.23 


5* 

1.23 


5* 

1.35 


5* 

1.20 


8 

1.29 


8 

1.27 


8 

1.28 


9* 

1.39 . 


9* 

1.37 


9* 

1.25 


15* 

1.31 


15* 

1.39 


15* 

1.24 


28 

1.43 


28 

1.38 


28 

1.20 


54 

■ 1.31 


52 : 

1.34 


52 

1.20 

5 

3 

1.29 

10 

3 

1.14 





3* 

1.29 


3* 

1.24 





5* 

1.30 


5* 

1.28 





8 

1.36 


8 

1.13 





9* 

1.38 


9* 

1.28 





15* 

1.48 


15* 

1.39 





28 

1.28 


26 

1.31 





54 

1.31 


52 

1.11 





*From bulk density data set obtained for the microwave remote. sensing 
experiment. 
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TABLE 2-2.- SOIL WATER PRESSURE VERSUS SOIL WATER CONTENT RELATIONSHIPS* 

FOR KEITH SILT LOAM 

[As reported by the Agricultural Technology Company] 

(a) Set I 


Field 

i 

Depth, 

in. 

Pressure 

i, -bars 

Field 

Depth, 

in. 

Pressure 

, -bars 

1/3 

15 

1/3 

15 

mnujuj, 

Water content, gg” 1 

i 

3 

wgm 

0.116 

8 

3 

mm 

0.116 


8 

mBSSm 



8 


0.105 


28 


KBe jH 


26 

0.265 

0.124 


54 

0.261 

0.127 


54 

0.292 


2 

3 

0.251 

0.108 

9 

3 

0.253 

0.106 


8 

0.259 

0.121 


8 

0.236 

0.011 


28 

0.262 

0.113 


28 

0.258 

0.122 


54 

0.263 

0.132 


52 

0.262 

0.115 

3 

3 

0.245 

0.104 

10 

3 

0.216 

0.089 


8 

0.236 

0.107 


8 

0.211 

0.091 


28 

0.258 

0.112 


26 

0.239 

0.102 


54 

0.262 

0.117 


52 

0.251 

0.115 

4 

3 

0.253 

0.112 

11 

3 

0.239 

0.098 

I 

8 

0.270 

0.133 


8 

0.247 

0.133 

i 

28 

0.277 

0.127 


28 

0.249 

0.128 

i 

54 

0.274 

0.122 


52 

0.254 

0.123 

! 5 

3 

0.267 

0.121 

12 

3 

0.270 

0.101 

i 

8 

0.257 

0.122 


8 

0.268 

0.112 . 

i 

28 

0.278 

0.124 


28 

0.259 

0.120 

i 

54 

0.278 

0.144 


52 

0.258 

0.118 

f 

6 

3 

0.281 

0.138 

13 

3 

0.262 

0.100 


8 

0.278 

0.141 


8 

0.259 

0.125 


24 

0.277 

0.144 


28 

0.257 

0.119 


54 

0.279 

0.153 


52 

0.262 

0.121 

7 

3 

0.240 

0.099 

14 

• 3 

0.280 

0.107 


8 

0.246 

0.118 


8 

0.262 

0.114 


28 

0.262 

0.122 


28 

0.268 

0.115 


52 

0.254 

0.116 


52 

0.269 

0.106 


♦Determined on crushed, sieved samples using a pressure plate apparatus. 
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TABLE 2-2.- Concluded, 
(b) Set II 


Field 

and 

location 

Depth, 

in. 


Pressure 

, -bars 


1/3 

D 

3 

6 

10 

15 


Water content, gg 

- 1 


2-3 

8 




0.160 

rm 

0.148 

2-3 

25 

0.328 

0.238 


0.160 


0.149 

2-3 

48 

0.279 

0.218 

0.161 

0.138 

B J 

0.119 

6-3 

8 

0.356 

0.274 

0.216 

0.203 

0.185 

0.183 

6-3 

25 

0.298 

0.217 


0.142 

0.139 

0.134 

6-3 

48 

0.284 

0.198 

0.141 

0.123 

0.117 

0.113 

'^11-4 

8 

0.328 

0.246 

0.198 

0.181 

0.139 

0.135 

11-4 

25 

0.326 

0.236 

0.187 

0.156 

0.151 

0.146 

11-4 

48 

0.279 

0.215 

0.158 

0.139 

0.124 

0.118 

14-1 

8 

0.298 

0.227 

0.176 

0.174 

0.147 

0.142 

14-1 

25 

0.309 

0.245 

0.193 

0.170 

0.135 

0.134 

14-1 

48 

0.279 

0.218 

0.159 

0.136 

0.130 

0.125 
















TABLE 2-3.- SATURATED HYDRAULIC CONDUCTIVITY (INCHES PER HOUR) 

IN KEITH SILT LOAM 

[As reported by the Agricultural Technology Company] 


Field 

and 

location 

I 

HH 

Flow time, hr 

Mean 


2 

4 

8 

24 

48 

All 

First 

three 

2-3 

8 

m 

E| 

0.46 

0.45 

0.45 

B 

0.41 

KB 

2-3 

25 

in 

EH 

1.08 

1.03 

1.07 

m 

1.07 

wm 

2-3 

48 

mm 

0.41 

0.48 

0.48 

0.46 

0.37 

0.47 

H 

6-3 

8 

B 

1.72 

2.06 

1.D4 

2.11 

1.51 

1.87 


6-3 

25 

EH 

0.22 

0.26 

0.26 

0.29 

0.33 

0.27 

m 

6-3 

48 

0.95 

0.79 

1.03 

1.00 

1.07 

1.08 

0.98 

0.92 

11-4 

8 

0.69 

0.55 

0.65 

0.55 

0.67 

0.77 

0.65 

0.63 

11-4 

25 

0.40 

— . 

0.43 

0.40 

0.46 

0.48 

0.43 

0.42 

11-4 

48 

0.43 

0.40 

0.52 

0.46 

0.48 

0.48 

0.46 

0.45 

14-1 

8 

0.72 

0.64 

0.77 

0.77 

1.03 

1.19 

0.85 

0.71 

14-1 

25 

1.38 

1.20 

1.46 

1.43 

1.44 

1.62 

1.42 

1.35 

14-1 

48 

0.41 

0.38 

0.43 

— 

— 

0.33 

0.39 

0.41 


Results: 

Overall mean = 1.90 cm/hr (0.748 in/hr) = 3.17xlO -2 cm/min 
Mean for A-horizon [20-cm (8-in.) depth] = 3.88xl0" 2 cm/min 
Mean for B-horizon [64-cm (25-in.) depth] = 3.25xl0 -2 cm/min 
Mean for C-horizon [122-cm (48-in.) depth] = 2.37xl0~ 2 cm/min 
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TABLE 2-5.- SOME IMPORTANT CHARACTERISTICS OF KEITH SILT LOAM, LOGAN COUNTY, KANSAS 



41-57 Cca -- - - — 0.1 11.4 65.7 22.8 53.7 23.5 0.21 1.20 












































Soil depth, in. 



Figure 2-1.- Bulk density profiles for Colby fields. Data for various 
fields are identified by field numbers. 
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Figure 2-2.- Particle-size distribution in Keith silt loam, Logan County, Kansas. Data represent 

the mean of two soil profiles. 




Organic carbon, percentage by weigh 



Figure 2-3.- Organic carbon distribution in Keith silt loam, Logan County, 
Kansas. Data from two soil profiles. The curve is hand drawn. 


2-12 




3. ESTIMATION OF HYDROLOGIC PROPERTIES 

3.1 SOIL WATER PRESSURE VERSUS SOIL WATER CONTENT RELATIONSHIPS 

Three different models were considered for estimating the relationship of soil 
water pressure to soil water content for Keith silt loam. These were: 

a. Regression model 

b. Rogowski model 

c. Ghosh model 

3.1.1 REGRESSION MODEL 

The regression model was developed by Hall et al. (ref. 9) in 1977 at the 
Rothamsted Experimental Station, England. A family of regression equations 
relates the volumetric water content at different pressures to particle-size 
composition, organic carbon, and bulk density. The equations are: 


Topsoils : 

e v (0.05) = 47.00 + 0.25C + 0.10Z + 1.12X - 16.52D b (1) 

e v (0.10) = 37.47 + 0.32C + 0.12Z + 1.15X - 12.5D b (2) 

e y (0.40) -=-26,66_ + 0.36C + 0.12Z + 1.00X - 7.64D b 

e v (2) = 8.70 + 0.45C + 0.11Z + 1.03X (4) 

6 v (15) = 2.94 + 0.83C - 0.0054C 2 (5) 

Subsoils : 

6 y (0.05) = 37.20 + 0.35C + 0.12Z - 11.73D b (6) 

6 V (0. 10) = 27.87 + 0.41C + 0.15Z - 8.32D b ' . (7) 

e y (0.40) = 20.81 + 0.45C + 0.13Z - 5.96D b (8) 

9 v (2) = 7.57 + 0.48C + 0.11Z (9) 

9 v (15) = 1.48 + 0.84C - 0.0054C 2 (10) 
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where 


e y = percentage of volumetric water content at different pressures (in bars) 
indicated in parentheses (One bar = 1000-centimeter height of a water 
column.) 

C = particles <0.002 millimeters, percentage by weight 
Z - particles 0.002 to 0.06 millimeters, percentage by weight 
X = organic carbon, percentage by weight 
D b = dry bulk density of the soil, grams per cubic centimeter 

For the various horizons of Keith silt loam, mean values of C, Z, X, and 
and applicable equations are as follows: 


Horizon 

C 

z i 

X 

Db 

Equations 

A 

28.5 

65.5 

1.30 

1.30 

1 through 5 

B 

31.0 

64.0 

— 

1.35 

6 through 10 

C 

22.0 

71.0 


1.28 

6 through 10 

Reference 

Fig. 2-2 

Fig. 2-2 

Fig. 2-3 

Fig. 2-1 



The relationships of soil water pressure to soil water content computed from 
equations (1) through (10) are plotted in figure 3-1. The soil water pres- 
sures for which water contents can be calculated from the regression equations 
are limited to -50 centimeters in the wet range and -15 000 centimeters in the 
dry range. Most soil moisture models and hydraulic conductivity calculations, 
however, require data beyond these limits. Thus, the utility of the regres- 
sion model appears to be limited. With this model, therefore, extrapolation 
of computed soil water pressure versus soil water content relationships may 
become necessary. In the dry range, the pressure versus water content rela- 
tionships for most soils, when plotted on semi log graphs, are very nearly 
linear, and extrapolation may not be too critical (see model comparisons by 
Rogowski, ref. 3; Baver et a!., ref. 10, p. 296). In the wet range, however, 
extrapolation may not be simple; and, as an aid, knowledge of other hydrologic 
parameters such as pressure and water content at the air-entry point and 
effective saturation may be required. For the purpose of comparative testing 
of soil moisture profile models, computed data were extrapolated to a water 
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content value equivalent to 85 percent of the theoretical porosity, 1 which 
appeared to approximate effective saturation. This was done merely to fill in 
the gap and is not recommended as a rule. 

The data in figure 3-1 show that the water retention characteristics of A and 
B horizons are similar and can be grouped together. Those for the C-horizon 
show substantially lower water retention. Thus, it appears that two distinctly 
different soil water pressure versus soil water content relationships, one for 
the combined A and B horizons and the other for the C-horizon, are applicable 
to Keith silt loam at Colby. This pattern of separation between horizons is 
also shown by the experimental data plotted in figure 3-1. 

Because experimental data were obtained on disturbed samples, only the 
pressure versus water content values in the dry range should be considered 
valid. In this range, a comparison between experimental and predicted data 
shows that, on an average, predicted water contents are 15 percent higher. A 
similar comparison of data in the wet range is not possible because (1) the 
regression model does not provide for computation of water contents in the wet 
range and (2) valid experimental data are not available. 

for ease in applying data to soil moisture models and in computing hydraulic 
conductivity, pressure-water content values for small increments of water 
content are presented in table 3-1. 

3.1.2 ROGOWSKI MODEL 

The Rogowski model (ref. 3) describes an analytic form of the soil water pres- 
sure versus soil water content relationship. The essential soil hydrologic 
parameters are water content at a soil water pressure of -15 000 centimeters, 
soil water pressure and water content at the air-entry point, and saturated 
water content. The model is of the following form. 

^f 8 [1 - (D b /D p )], where f is the theoretical porosity, D b is the bulk 

density, and 0 p is the particle density assumed to be 2.65 grams per cubic 
centimeter. y 
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e = (0 1 v a An(i|> * + 1 ) ■ ; ip 2 > ij> 
v ' v ' e e . e 


e v = < 0 v>e 


♦ < ^ 


( 11 ) 

( 12 ) 


° = C( 0 v )i5 “ ( 0 v ) e ]/^n(i{» i5 - ip e + 1) (13) 

v/here e v is the volumetric water content in cubic centimeter per cubic centi- 
meter and i|> is the soil water pressure in centimeters (in terms of magnitude 
only). Subscripts e and 15 denote air entry and 15 000 centimeters pressure, 
respectively. 


For Keith silt loam at Colby, experimentally measured ( e v )^5 values (i.e., 
volumetric water content at a pressure of -15 000 centimeters) can be 
accepted. Values are given in table 2-2 and are plotted in figure 3-1. 
Parameters (8 ) and ij> e were not measured and must be estimated based on data 
in the literature. Some data compiled by Rogowski (ref. 11) are summarized in 
table 3-2. These data show that, on an average, the ratio of air-entry water 
content to saturated water content is approximately 0.9 and the air-entry 
pressure is approximately 23 centimeters. Assuming these values are appli- 
cable to the Keith silt loam at Colby, relevant hydrologic parameters for 
computing pressure-water content relationships can be developed. Important 
parameters for Keith silt loam are shown below. 


Horizon 

V 

gem - 3 

Db> 

gcm - 3 


(®v)e/ ( 8 9 v)q 

(9y)e> 

3 -3 

enrem J 

V 

cm 

( 9 v ) 1 5 » 

3 -3 

enrem J 

-a 

A 

2.65 


0-.509 

0.90 

0.453 

23.0 

0.198 

0.02704 

B 

2.65 

■ 

0.491 

0.90 

0.442 

23.0 


0.02621 

C 

2.65 


0.517 

0.90 

0.465 

23.0 

0.152 

0.03255 


Volumetric water contents as a function of soil water pressure, computed from 
equations (11) through (13), are shown in figure 3-2. These equations allow 
computation of water content at pressures below. the air-entry point only. For 
soil water pressures above the air-entry point, Rogowski suggests an alterna- 
tive formulation in the form: 

8 v = ( 9 v^e + 8 Zn (^ e " * + *) l * * + e ( 14 ) 
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«, > <\) e US) 

6 - C(e v ) 0 * (»,) e )/«»(* e -*„ + x) (16) 

where 4» Q and (0 y ) are saturation soil water pressure and water content, 
respectively. By using equations (14) through (16), the pressure-water 
content relationships were extended to the saturation point. 

The Colby experimental data are shown also in figure 3-2. Agreement between 
experimental and predicted water contents appears reasonable in the pressure 
range of -3000 to -15 000 centimeters. At higher pressures, experimental data 
show significantly higher water contents. However, because Colby data were 
obtained on disturbed samples, it would be inappropriate to use them for 
verifying the predicted results. Rogowski compared the predicted and experi- 
mental soil water pressure versus soil water content relationships for a 
number of soils and concluded that, in view of field variabilities, results 
obtained from equations (11) through (16) are quite adequate. 

Table 3-3 shows the pressure-water ,content values for small increments of 
water content. 

3.1.3 GHOSH MODEL 

The Ghosh model (ref. 12) describes the pressure-water content relationship of 
soils in the following manner: 

♦ - v e [o v /(9 v ) 0 r s - - . ■ (l?) 

where > is any soil water pressure, 4> e is the soil water pressure at air 
entry, e v is the water content corresponding to soil water pressure ( 0 V ) O 
is the water content at saturation, and 6 is an empirically determined 
constant. For application of this model, at least one measurement of 6 y at ^ 
and knowledge of i|> e , particle-size distribution, and bulk density are 
necessary. The value of (8 V ) 0 can be estimated from a knowledge of bulk 
density, while B is given by the following formulation: 

B = 26.5(X 2 /X 1 ) 1 * 786 (18) 
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where X 1 and X 2 are the percentages by weight of sand 2 and silt 3 content of 
the soil. 

For Keith silt loam, bulk density and particle-size distribution can be 
inferred from figures 2-1 and 2-2, respectively. 

Measured water content values at pressures ranging from -3000 to -15 000 centi- 
meters are accepted as being reliable. These values can be inferred from 
figure 3-1 or 3-2. 


Relative to the Ghosh model, various parameters for Keith silt loam are as 
follows: 


Horizon 

V 

_3 

gem J 

Db» 

_3 

gem ° 

(8 V ) 0 . 

cnrcnT^ 

(9 v )l5» 

3 -3 

cm J cm 0 

* 4 
v e* 

-cm 

X 1 » 
% 

x 2* 

% 

8 5 

A 

2.65 

1.30 

0.509 

0.198 

255 

52.5 

19.0 

4.314 

B 

2.65 

1.35 

0.491 

0.190 

94 

49.0 

20.0 

5.348 

C 

2.65 

1.28 

0.517 

0.152 

11 

54.5 

23.5 

5.899 


The computed values of for the A and B horizons, although not impossible, 
appear too high when compared with published experimental data (see table 3-2). 
The empirical formulation for 8 in the Ghosh model [equation (18)] was 
developed from data on sandy materials; it may not be adequate for heavier 
textured soils. Furthermore, variabilities associated with soil density, 
estimates of saturated water content, and the choice of a pair of measured 
4»-0 v values as inputs to the model should be expected to have an influence on 

v , .. 

2 ; 

Particles ranging in size from 2.0 to 0.02 millimeters. 

3 

Particles ranging in size from 0.02 to 0.002 millimeters. 

4 

Computed on the basis of water content at a pressure of -15 000 centimeters, 

8, and (0 V ) O as shown. 

Computed on the basis of particle-size parameters and x 2 » 
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An alternative approach is to estimate both i|> and 0 from a pair of 
simultaneous equations. However, this approach requires at least two pairs of 
measured 4»-e v values, in addition to saturated water content. 

I 

For Keith silt loam at Colby, measured water contents at soil water pressures 
of -3000 and -15 000 centimeters (see figs. 3-1 and 3-2) were used for esti- 
mating tg and 0 from equation (17). The substitution of appropriate values 
for 0 v » and ( 0 V ) O in equation (17) resulted in the following revised values 

of and e for the various horizons, 
e 


Horizon 

V 

-cm 

0 

A 

37 

6.362 

B 

44 

6.134 

C 

9 

6.088 


Computed soil water pressure versus soil water content relationships based on 
revised i|> e and 0 values and equation (17) are shown in figure 3-3. The data 
in figure 3-3 show that the A and B horizons have very nearly identical soil 
water pressure versus soil water content relationships and that water retention 
for these two horizons is substantially higher than 'for-ihe C-horizon. 

Table 3-4 shows the pressure-water content values for small increments of 
water content . 

3.2 HYDRAULIC CONDUCTIVITY VERSUS SOIL WATER CONTENT RELATIONSHIPS 

A knowledge of hydraulic conductivity as a function of soil water content is 
important in analyses of water movement. Although experimental data obtained 
in the field or on undisturbed core samples are preferred, soil variability is 
often so large that experimental efforts needed for large-scale application 
become prohibitive. As a result, predictive methods have gained popularity; 
in many instances (e.g., Nielsen et al., ref. 4), satisfactory estimates of 
hydraulic conductivity have been obtained. 


3-7 




For Keith silt loam at Colby, saturated hydraulic conductivity for selected 
depths in four fields was measured. These data are reported in table 2-3. 
Unsaturated hydraulic conductivity was not measured but was predicted by 
Jackson's (ref. 13) modification of Marshall's (ref. 14) pore-interaction 
model. 6 The necessary input for Jackson's method is the soil water pressure 
versus soil water content relationship for the soil. The curve describing 
this relationship is divided into n equal water content increments, and the 
conductivity is calculated by the following equation: 




n r 

L, (2j + 1 

- 2i)*- 2 

j=i L 

J 

n r 


r (2j - 


j=i L 

J J 



(19) 


where 

K.j = hydraulic conductivity (centimeters per minute) of water content at the 
upper end of the i^ water content interval 

K $ - measured saturated hydraulic conductivity (centimeters per minute) 

e.^water content at the upper end of the i th water content interval (cubic 
centimeter per cubic centimeter) 


= highest or saturated water content (cubic centimeter per cubic 
centimeter) 

p = an empirical constant [The value depends on the method of computation; 
for equation (19), p = 1.] 

ij/j = soil water pressure at the midpoint of water content interval (minus 
centimeters of water) 


The computer program for computation of unsaturated hydraulic conductivity is 
given in the appendix. 
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For the combined A and B horizons of Keith silt loam, hydraulic conductivity 
data computed from equation (19) are shown in figure 3-4. Those for the 
C-horizon are shown in figure 3-5. Hydraulic conductivities for each of the 
two sections of the soil profile are based on soil water pressure versus water 
content relationships predicted from the regression, Rogowski, and Ghosh 
models. Table 3-5 shows the hydraulic conductivities for small increments of 
water content. 
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TABLE 3-1.- SOIL. WATER, PRESSURE VERSUS SOIL WATER CONTENT .RELATIONSHIPS* 
FOR KEITH SILT LOAM COMPUTED FROM THE REGRESSION MODEL 


" ’ Soil' horizon ' j 

A + B [0 to 102-cm (40-in.) depth] 

C [bel ow 102-cm 

(40-in.) depth] j 

Water 

content, 

3 -3 

ctrrcm ° 

Pressure, 

-cm 

Water 

content, 

cnr^cm -3 

• Pressure* 
-cm 

• 

, Water 
content, 

cnPcm" 3 

Pressure, 

-cm 

• 

Water 

content, 

cn^cm -3 

p I 

Pressure, 

-cm 

0.43 

1.00x10° 

0.23 

1.30xl0 4 

0.43 

1.00x10° 

0.23 

4.20x10° | 

0.42 

9.00x10° 

0.22 

1.65xl0 4 

0.42 

4.00x10° 

0.22 

5.40xl0 3 | 

0.41 

3-OOxlO 1 

0.21 

2.25xl0 4 

0.41 

1.30X10 1 

0.21 

6.80xl0 3 ! 

0.40 

6.60X10 1 

0.20 

2.95xl0 4 

0.40 

2.50X10 1 

0.20 

8.60x10° 

0.39 

1.20xl0 2 • 

0.T9 

•3.80xl0 4 < 

0.39 

4.20X10 1 

0.19 

1.07xl0 4 

0.38 

1.90xl0 2 

0.18 

4.90xl0 4 

0.38 

6.30X10 1 . 

0.18 

1.35xl0 4 

0.37 

2.60xl0 2 

0.17 

6.40xl0 4 

0.37 

9.60X10 1 

0.17 

1.70xl0 4 

! 0.36 

3.60xl0 2 

0.16 

8.60xl0 4 

0.36 

1.35xl0 2 

0.16 

2.10xl0 4 

0.35 

5.00xl0 2 

0.15 

1.15xl0 5 

0.35 

1.90xl0 2 

0.15 

2.60xl0 4 

0.34 

6.60xl0 2 

0.14 

1.50xl0 5 

0.34 

2.55xl0 2 

0.14 

3. 30x1 0 4 

0.33 

7.70xl0 2 

0.13 

1.95xl0 5 

0.33 

3.30xl0 2 

0.13 

4.10xl0 4 

0.32 

1.15x10° 

0.12 

2.50xl0 5 

0.32 

4.35xl0 2 

0.12 

5.00xl0 4 ! 

0.31 

1.50xl0 3 

0.11 

3.40xl0 5 

0.31 

5.60xl0 2 

0.11 

6.40xl0 4 

0.30 

2.00xl0 3 

0.10 

4.40xl0 5 

0.30 

7.40xl0 2 

0.10 

8.00xl0 4 

0.29 

2.50xl0 3 

0.09 

5.70xl0 5 

0.29 

9.60xl0 2 

0.09 

9.80xl0 4 

0.28 

3.30xl0 3 

0.08 

7.40xl0 5 

0.28 

1.23xl0 3 

0.08 

1.23x10° 

0.27 

4.30xl0 3 

0.07 

9.50xl0 5 

0.27 

1.60xi0 3 

0.07 

1.55x10° 

0.26 

5.70xl0 3 

0.06 

1.30x10° 

0.26 

2.00xi0 3 

0.06 

1.95x10° 

0.25 

7.50xl0 3 

0.05 

I.65xl0 6 

0.25 

2.60xl0 3 

0.05 

2.45x10° 

0.24 

9.70xl0 3 

0.04 

2.20xl0 6 

0.24 

3.30xl0 3 

- - — 

0.04 

3.00x10° 
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TABLE 3-2.- NUMERICAL VALUES OF IMPORTANT PHYSICAL AND 
HYDROLOGIC PARAMETERS FOR SEVERAL SOILS 

[From prepublication copy of Watershed Physics: Moisture 
Characteristics and Variability Criteria by 
A. S. Rogowski , USDA Agricultural Research Service, 
Beltsville, Maryland] 

(a) 15-bar* and air-entry values 


Soil 

series 

^VlS* 

cnAcm"^ 

V* 

-cm 

( 9 v ) e . 

crti^cm*^ 

Watson S. 

— 

39.8 

0.350 

Adelanto c.l . 

0.174 

12.0 

0.414 

Adelanto c .1 . 

0.110 

10.0 

0.402 

.Panoche c.l. 

0.213 

35.0 

0.448 

Panoche c .1 . 

0.275 

25.0 

0.431 

Panoche c.l . 

0.199 

32.0 

0.428 

Panoche c .1 . 

0.196 

47.0 

0.476 

. Panoche c.l . 

— 

32.0 

0.491 

Panoche c .1 . 

— 

32.0 

0.522 

Yolo 1. 

0.167 

10.0 

0.482 

Yolo 1. 

0.167 

20.0 

0.482 

Miller SI .c. 

— 

10.0 

0.357 

Miller Si .c. 

0.200 

15.0 

0.340 

Miller Si .c. 

0.135 

30.0 

0.333 

Miller Si.c. 

0.112 

30.0 

0.340 

Monona Si .1 

0.125 

9.0 

0.476 

Adams F .s .1 . 

0.102 

10.0 

0.519 

Cecil C. 

0.299 

100.0 

0.422 

Houston C. 

0.317 

200.0 

0.428 


(b) Density and C(0 v ) e /(9 y ) o 3 ratios 


Soil 

series 

V 

gem' ^ 

K3| 

1 

( 9 v)o* 

cirAcnf 

. (®y)o» 
cnPcnf ^ 

C(e v ) e /(9 v ) 0 D ? 

Panoche 

2.65 

1.31 

mm 

mm 

mgm 

Panoche 

2.65 

1 .25 

■El 

IBS 


Adel anto 

2.72 

1.47 

0.460 

|E9 

0.89 

Yolo 

2.65 

1.23 

0.536 

ESI 

0.89 

Mi 1 1 er 

2.65 

1.61 

0.393 

KE9 

0.85 

Webster 

2.58 

1.28 

0.504 

■IS 

• 0.83 

Webster 

2.58 

1.20 

0.535 

0.513 

0.96 


•Represents a soil water pressure of -15 000 on. 
*Mean = 23 cm for all soils except sand and clay. 
Iflean C(9 V ) e /(9 v ) 0 3 = 0.90. 
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TABLE 3-3.- SOIL WATER PRESSURE VERSUS SOIL WATER CONTENT RELATIONSHIPS* 
FOR KEITH SILT LOAM COMPUTED FROM THE ROGOWSKI MODEL 


Soil horizon 


■ .. 

A + B [0 to 102-cm (40-in.) depth] 

C [below 102-cm 

(40-in.) depth] 

Water 

content, 

i cm 3 cnf 3 

Pressure, 

-cm 

Water 

content, 

cnr^cm" 3 

Pressure, 
• -cm 

Water 

content, 

cm 3 cm" 3 

Pressure, 

-cm 

Water 

content, 

cnPcm" 3 

Pressure, 

-cm 

1 0.50 

1.00x10° 

0.27 


0.52 

1.00x10° 

0.27 

4.30xl0 2 

0.49 

1.20X10 1 

0.26 

1.35xl0 3 

0.51 

9.00x10° 

0.26 

5.80xl0 2 

0.48 

1.80X10 1 

0.25 

2.00xl0 3 * 

0.50 

1.60X10 1 

0.25 

7.90xl0 2 

0.47 

2.00x10* 

0.24 

2.85xl0 3 

0.49 

.1.95*10* 

0.24 

1.03xl0 3 

0.46 

2.20X10 1 

0.23 

4.15xl0 3 

0.48 

2.10x10* 

0.23 

1.45xl0 3 

0.45 

2.30xl0 x 

0.22 

6.00xl0 3 

0.47 

2.20X10 1 

0.22 

1.90xl0 3 

0.44 

2.40X10 1 

0.21 

8.60xl0 3 

0.46 

2.30X10 1 

0.21 

2.60xl0 3 

0.43 

2.60X10 1 

0.20 

1.25xl0 4 

0.45 

2.40X10 1 

0.20 

3.60xl0 3 

0.42 

2.80X10 1 

0.19 

1.80xl0 4 

0.44 

2.53X10 1 

0.19 

4.80xl0 3 

0.41 

3-OOxlO 1 

0.18 

2.60xl0 4 

0.43 

2.65X10 1 

0.18 

6.50xl0 3 

0.40 

3.30X10 1 

0.17 

3.80xl0 4 

0.42 

2.80x10* 

0.17 

8.90xl0 3 

0.39 

3.60X10 1 

0.16 

5.37xl0 4 

0.41 ' 

3.00x10* 

0.16 

1.20xl0 4 

0.38 


0.15 

7.81xl0 4 

0.40 


0.15 

1.65xl0 4 

0.37 

4.60X10 1 

0.14 

1.14x105 

0.39 

3.50x10* 

0.14 

2.25xl0 4 

0.36 

5.30X10 1 

0.13 

1.65x10° 

0.38 

3.85x10* 

0.13 

3.00xl0 4 

0.35 

6.50X10 1 

0.12 

2.41x10° 

0.37 

4.30x10* 

0.12 

4.10xl0 4 

0.34 

8.60X10 1 

0.11 

3.51xl0 5 ' 

0.36 


0.11 

5.40xl0 4 

0.33 

1.15xl0 2 

0.10 

5.11x10° 

0.35 



7.41xl0 4 

0.32 

1.55xl0 2 

0.09 

7.43xl0 5 

0.34 

6.80x10* 

0.09 

1.01x10° 

; 0.31 

2.20xl0 2 

0.08 

1.08x10° 

0.33 

8.60x10* 

0.08 

1.37x10° 

j 0.30 

3.15xl0 2 

0.07 

1.57x10° 

0.32 

l.lOxlO 2 

0.07 

1.86x10° 

! 0.29 

4.60xl0 2 

0.06 

1.57x10° 

0.31 

1.40xl0 2 

0.06 

2.53x10° 

0.28 

6.60xl0 2 

0.05 

3.34x10° 

0.30 

1.85xl0 2 

0.05 

3.44x10° 



0.04 

4.86x10° 

0.29 

2.40xl0 2 

0.04 

4.68x10° 





0.28 

3.20xl0 2 




♦Interpolated from figure 3-2. 
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TABLE 3-4.- SOIL WATER PRESSURE VERSUS SOIL WATER CONTENT RELATIONSHIPS* 
FOR KEITH SILT LOAM COMPUTED FROM GHOSH MODEL 


Soil horizon - • 

A + B [0 to 102-cm (40-in.) depth] 

C [below 102-cm (40-in. j depth] 

Water 

content, 

cm 3 cm“ 3 

Pressure, 

-cm 

Water 

content, 

cnr^cm" 3 

Pressure, 

-cm 

Water 

content, 

cnr^cm" 3 

Pressure, 

-cm 

Water: 

content, 

cnPcnf 3 

Pressure, 

-cm 

mm 

1.00x10° 

mm 

1.95xl0 3 . 

0.52 

1.00x10° 

0.27 

4.80xl0 2 

■n 

4.50X10 1 


2.45xl0 3 

0.51 

9.80x10° 

0.26 

6.00xl0 2 

0.48 

5.20X10 1 

■9 

3.20xl0 3 

0.50 

l.lOxlO 1 

0.25 

7.70xl0 2 

0.47 

6.00X10 1 

0.24 

4.00xl0 3 

0.49 

1.24X10 1 

0.24 

9.80xl0 2 

0.46 

6.80X10 1 

0.23 

5.30xl0 3 

0.48 

1.40X10 1 

0.23 

1.25xl0 3 

0.45 

7.90xl0 ! 

0.22 

6.80xl0 3 

0.47 

1.60X10 1 

0.22 

1.65xl0 3 

0.44 

9-lOxlO 1 

0.21 

9.00xl0 3 

0.46 

1.84xl0 x 

0.21 

2.20xl0 3 

0.43 

1.05xl0 2 

0.20 

1.25xl0 4 

0.45 

2.10X10 1 

0.20 

3.00xl0 3 

0.42 

1.22xl0 2 

0.19 

1.74xl0 4 

0.44 

2.40X10 1 

0.19 

4.10xl0 3 

0.41 

1.40xl0 2 

0.18 

2.40xl0 4 . 

0.43 

2.75xlO x 

0.18 

5.80xl0 3 

0.40 

1.64xl0 2 

0.17 

3.60xl0 4 

0.42 

3.15xlO x 

0.17 

8.20xl0 3 

0.39 

1.92xl0 2 

0.16 

5.30xl0 4 

0.41 

3.65X10 1 

0.16 

1.15xl0 4 

0.38 

2.25xl0 2 

0.15 

7.58xl0 4 

0.40 

4.30X10 1 

0.15 

1.65xl0 4 

0.37 

2.70xl0 2 

0.14 

1.17xl0 5 

0.39 

5.00xl0 x 

0.14 

2.55xl0 4 

0.36 

3.20xl0 2 

0.13 

1.85xl0 5 

0.38 

5.80xl0 x 

0.13 

4.15xl0 4 

0.35 

3.80xl0 2 

0.12 

3.06xl0 5 

0.37 

T.OOxlO 1 

0.12 

6.54xl0 4 

0.34 

4.60xl0 2 

0.11 

5.27xl0 5 

0.36 

8-OOxlO 1 

0.11 

l.llxlO 5 

0.33 

5.50xl0 2 

0.10 

9.55xl0 5 

0.35 

9.70xl0 x 

0.10 

1.99x10® 

0.32 

6.70xl0 2 

0.09 

1.84x10® 

0.34 

1.15xl0 2 

0.09: 

3.77x10® 

0.31 

7.40xl0 2 

0.08 

3.85x10® 

0.33 

1.40xl0 2 

0.08 ; 

7.73x10® 

0.30 

l.OOxlO 3 

0.07 

8.87x10® 

0.32 

1.68xl0 2 

0.07 

1.74x10® 

0.29 

1.25xl0 3 

0.06 

2.32xl0 7 

0.31 

2.04xl0 2 

: 0.06 

4.45x10® 

0.28 

1.55xl0 3 

0.05 

7.26xl0 7 

0.30 

2.50xl0 2 

; o.p5: 

1.35xl0 7 



0.04 

2.93xl0 8 

0.29 

3.05xl0 2 

0.04 

5.26xl0 7 





0.28 

3.85xl0 2 . 


• 


interpolated from figure 3-3. 
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TABLE 3-5.- SOIL WATER CONTENT VERSUS HYDRAULIC CONDUCTIVITY* RELATIONSHIPS FOR KEITH SILT LOAM 
BASED ON DIFFERENT MODELS FOR PRESSURE-WATER CONTENT RELATIONSHIPS 
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Interpolated from figures 3-4 and 3-5. 



Soil water pressure, -cm H 2 0 



Figure 3-1.- Soil water pressure versus soil water content relationships 
for Keith silt loam computed from the regression model . 
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4. CONCLUDING REMARKS 


For hydraulic conductivity versus soil water content relationships, Marshall's 
model and its modifications have been shown to represent field conditions 
adequately. However, because the soil water pressure versus water content 
relationship and saturated hydraulic conductivity are key inputs to the 
hydraulic conductivity model, the accuracy of these properties would largely 
determine the quality of the predicted hydraulic conductivity data. 

The regression, Rogowski, and Ghosh models, used to obtain soil water pressure 
versus soil water content relationships for Keith silt loam at Colby, repre- 
sent two approaches. In the regression- type approach, no attempt is made to 
incorporate hydrophysical characteristics of the soil into the model. Hence, 
predicted results may be useful only to evaluate variations between soils 
resulting from such factors as texture, bulk density, and organic matter. The 
Rogowski and Ghosh models, on the other hand, represent the analytical 
approach and do, in fact, consider some well -recognized hydrophysical quan- 
tities, such as water content at one or two specified soil water pressures, 
soil water content and pressure at the air-entry point, and saturated water 
content, thus, the results of the analytical approach may represent a more 
realistic description of soil hydrologic characteristics than wi4T-the_ regres- 
sion approach. Colby experimental water contents in the drier range were used 
as inputs to the Rogowski and Ghosh models; therefore, agreement between 
experimental and predicted data in that range does not constitute proof of the 
models' predictive ability. Additional tests over a wide range of water con- 
tents need to be carried out. In the wet range, soil water pressure and water 
content corresponding to the air-entry point and saturation values were based 
on data in the literature and may not be representative of Keith silt loam 
soils at Colby. Thus, the reliability of data generated from the Rogowski and 
Ghosh models should be considered within the limits of the assumptions that 
experimental water contents in the drier range are valid and that the 
literature-based air-entry and saturation values are applicable to Keith silt 
loam. 
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5. RECOMMENDATIONS 

Analyses of physical and hydrologic properties of Keith silt loam indicate 
similarities between the A and B horizons, which together appear distinctly 
different from the C-horizon. As far as possible, the two sections of the 
profile should be treated separately. 

* ” '• . • • . 

I . 

Computed hydrologic properties presented in this report appear typical of most 
loam to silt loam soils and should be adequate for the- purpose of comparing 
different soil moisture models. The results from analytical models should be 
preferred over those from regression-type models. 

Concerning the effects of hydrologic properties on the ability of a soil 
moisture model to predict the measured soil moisture regime, experimentally 
determined hydrologic properties would appear more appropriate than those 
predicted. When detailed information is not available, experimental determi- 
nation of the soil water pressure versus soil water content relationship near 
saturation, at the air-entry point, and at a pressure of -15 000 centimeters 
should be the minimum requirement. In the wet range, field measurements 
should be preferred over laboratory measurements. 

In the application of soil moisture models to the task of predicting soil 
moisture regime over large areas comprising a variety of soils, it is highly 
unlikely that experimental data pertinent to different soils would be avail- 
able. Nor would it be economically feasible to devote experimental efforts. 
Furthermore, in view of the large variability between soils (even within a 
series), such efforts would not be justified. Therefore, much of the applica- 
tion of soil moisture models will have to depend on predicted hydrologic 
properties. 
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In order to develop reasonable predictive methods, it is recommended that 

a. Models for predicting soil hydrologic properties be extracted from the 
1 iterature 

b. For as many soils as possible, a catalog of experimentally determined soil 
hydrologic properties and such basic properties as are necessary as inputs 
for the models be developed, maintained, and updated 

c. The effects of predicted soil hydrologic properties on predictive 
performance of soil moisture models be evaluated 

d. Efforts be directed toward modification and improvement of existing models 
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APPENDIX 

FORTRAN PROGRAM FOR THE COMPUTATION OF 
UNSATURATED HYDRAULIC CONDUCTIVITY 


FILE : HYCOND F OH TRAN A 


CUNVtPSAT IONAL MON I TOW SYSTEM 


HYCONO PPOOpAM 

THIS PpUOpAM CALCULATES HYDRAULIC CONDUCTIVITY FROM SOILS DATA. 
EPSlLN IS WAThP CUNlENT IN VOLUME I Wit FRACTION. 

SKAY IS MEASuplU SAT UP 4 f EU nYUHAULlC CUNUUC T 1 V 1 T Y . 

M IS THE NUMhEh oh ..ATtH CuwIEnI INTERVALS POM WHICH 
COMPUTATIONS APt. MADE. 

H IS Tut SU1L v* A T t P ppESSUhE AT THE MIUHOINT OF THE WATER 
CONTENT INTERVALS. 

TaaY IS T fit CALCULATED HYOHAULlC CONDUCTIVITY. 

DIMENSION n(o(j)« EPSiLNlSo) 

REA0<S»2S> SPAY.P.M 
FOHMAT (2A»Ell.t»FS.2»U) 

WHITE (to«2b) siv a Y » P » M 
DO 5 1=1 . m 

REAUtb. 30) H(I) .EPsILNCI) 

WHITE (6*30) H(l) .EPSILnUJ 
FORMAT (lA.El2.AtFA. 2) 

SUM2=0 

OU 20 1=1 . m 

SUMl=u 

DO 10 J= I ♦ M 

X=I 

Y “ J 

VM.U£1=<2*Y*1-2*X)/H(J>**2 
IF ( I .OT. I ) 00 TU 7 
VALUE2=<2*Y-1)/h(J)**2 

SUM2=SUM2* V ALUE2 

Sum i =sumi .value 1 

wkITE ( 16. 3o) value 1 .value 2 

CONTINUE 

FAC TH=SUM1/SUm2 

TKAY=(SKAYJ*( (EPSILN<l)/£PSlLNU))**P)*FACTR 

wRITE(6«4i) SUM i . SUM2. Tin AY 

PUNCH4 I * SUM I » SUM2 . Tf\ a Y 

WK IT E ( I to » 4 1 ) SUM 1 . SUM2. TNAY 

FOpmaT ( 1A. JE12.4./) 

CONTINUE 

STOP 

ENU 
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